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Intercomparisons  of  Radiosondes  and  an  Airborne  Refractometer 

for  Measuring  Radio  Ducts 


1.  INTKODUCTION 

In  September  1985.  a  series  of  field  tests  was  undertaken  at  Chatham,  MA  as 
part  of  an  in  house  program  to  study  atmospheric  effects  on  Air  Force  microwave 
communication  systems.  Part  of  this  study  used  refractive  index  (N)  profile  data 
from  three  sources:  an  aircraft  refractometer,  a  portable  radiosonde,  and  a 

synoptic  NWS  rawinsonde.  The  primary  purpose  of  this  series  of  tests  was  to 

2  ) 

measure  the  refractive  index  structure  function  (C^  in  and  near  atmospheric 

inversion  layers.  The  secondary  purpose  was  to  compare  refractive  index  profile 

data  from  all  three  sources.  This  report  presents  only  those  data  pertaining  to  the 

secondary  purpose.  Some  of  the  shortcomings  of  the  data  base,  such  as  incomplete 

2 

aircraft  profiles,  are  due  to  the  precedence  of  the  measurements  in  flight 
planning. 

Chatham  is  located  on  the  southeastern  tip  of  Cape  Cod.  It  was  chosen  because 
of  the  presence  of  a  National  Weather  Service  upper  air  station,  its  proximity  to  our 
work  location,  its  proximity  to  the  ocean,  and  the  presence  of  an  airport. 


(Keeeived  for  publication  1  July  1906) 


2.  DESCRIPTION 


Three  methods  were  used  to  obtain  vertical  profiles  of  the  microwave  refrac' 
tive  index.  Two  of  these  are  radiosonde  systems  that  measure  temperature,  pres¬ 
sure,  and  humidity,  and  the  third  is  a  refractometer  that  measures  index  of 
refraction  directly. 

The  first  radiosonde  system  is  called  the  Atmospheric  Data  Acquisition  System 
(ADAS).  It  consists  of  a  ground  receiver  and  data  processor,  a  cassette  tape 
recorder,  a  printer,  and  the  radiosonde,  as  shown  in  Figure  1.  The  radiosonde 
itself  consists  of  a  bead  thermistor  to  measure  the  temperature,  a  capacitively 
sensed  aneroid  cell  to  measure  the  pressure,  and  a  carbon  hygristor  to  measure 
the  humidity.  The  sonde  time-commutates  the  data  from  these  sensors  and  uses 
it  to  modulate  a  crystal  controlled  transmitter  at  403.5  MHz.  The  commutation 
provides  a  complete  frame  (T,  P,  H)  about  every  5  sec.  All  of  these  data  were 
recorded/printed  and  were  available  for  analysis. 


Figure  1.  Portable  Radiosonde  System  -  Automatic  Data  Acquisition  System 
(ADAS)  -  Manufactured  by  AIR  Inc.  ,  Boulder,  CO 
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The  second  radiosonde  system  is  used  at  synoptic  upper  air  stations  in  the 
U.  S.  *  The  radiosonde  (VIZ  mfg.  No.  1392-510)  consists  of  a  rod  thermistor  to 
measure  temperature,  an  aneroid  cell  to  measure  pressure,  and  the  same  carbon 
hygristor  to  measure  humidity.  This  radiosonde  commutates  the  temperature  and 
humidity  data  along  with  references  and  send  them  to  the  ground  station  as  modula¬ 
tion  on  a  1680  MHz  carrier  frequency.  A  principal  difference  between  the  sondes 
is  that  the  aneroid  cell  drives  a  mechanical  commutator  in  the  synoptic  instrument. 
Consequently,  the  position  of  the  commutator,  or  contact  number,  contains  the 
pressure  information.  The  commutation  rate  for  the  synoptic  sonde  is  thus  deter¬ 
mined  by  the  rise  rate,  and  is  slower  than  the  commutation  rate  for  the  ADAS 
sonde.  The  time  for  each  contact  averages  about  20  sec  in  the  lower  atmosphere 
and  the  cycle  is  THTHTHTHTHTHTHTIITHTKT  ...  where  T  is  temperature,  11  is 
relative  humidity,  and  K  is  reference.  This  results  in  a  time  between  adjacent 
humidity  values  of  20  to  30  sec,  and  when  a  reference  intervenes,  of  about  50  sec. 

The  refractometer  used  was  developed  at  the  National  Telecommunications 

2 

and  Information  Sciences  (NTIS)  Laboratory  at  Boulder,  CO.  It  is  a  single  cavity 
Birnbaum  type.  The  cavity  is  not  the  frequency  determining  element,  but  is 
excited  with  a  signal,  centered  on  the  resonant  frequency,  that  has  a  sawtooth 
frequency  modulation.  An  error  signal  that  detects  when  the  resonance  is  not 
centered  is  used  to  adjust  the  excitation.  The  refractometer  relies  on  the  relation¬ 
ship 


AN 


Af 

f 


10 


(1) 


so  that  the  resonant  frequency  (f)  of  the  cavity  depends  on  the  refractive  index 

(N)  of  the  air  that  is  continuously  flushed  through  it.  The  refractometer  is  mounted 

on  the  str  ut  below  the  wing  of  a  single  engine  aircraft,  as  shown  in  Figure  2.  The 

refractive  index,  together-  with  the  air  temperature  from  a  bead  thermistor  and 

atmospheric  pressure  from  u  laboratory  aneroid  cell,  are  recorded  twice  a  second 

on  magnetic  tape.  Although  recorded  twice  a  second,  the  refractive  index  from  the 

refractometer  is  updated  only  once  a  second.  The  aircraft  recording  is  controlled 

by  a  computer  system  that  digitizes  the  data  (12  bit)  prior  to  recording  it,  and  also 

allows  the  operator  to  insert  comments  on  atmospheric  conditions  and  the  tvpe  of 
2 

test.  I  or  the  C ^  measurements,  additional  data  are  recorded  at  high  speed 
(6-1  samples/sec)  from  a  second  channel  of  the  refractometer. 


1 .  Federal  Meteorological  Handbook  No.  3  (1971)  Radiosonde  Observations  - 

Supe  rintendent  of  Documents,  U.  S.  Gov't  Printing  (Office,  Washington,  D.f. 

2.  Thompson,  M.  C.  ,  Marler,  F.  E.  ,  and  Allen,  K.  I..  (  1980)  Measurement  of  the 

microwave  structure  constant  profile,  IEEE  Trans.  Antennas  I’ropag.  , 
AP-28(No.  2).  '  ^ 
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1  igure  2.  Aircraft  Kefractometer  Configuration  -  Insert  Shows  Internal 
Kleetronic  Equipment 


:t.  FISTS 


All  radiosonde  tests  wore  performed  at  the  NWS  upper  air  station  at  Chatham, 

M  \.  Tile  station  is  located  on  a  bluff  1 1>  meters  above  the  ocean.  The  re  were  two 
tapes  of  soundings,  those  involving  both  the  synoptic  radiosonde  and  the  AI>AS  sonde, 
in  1  those  involving  only  the  ADAS  sonde.  In  the  double  sounding  case  the  ADAS 
sonde  was  piggy-hacked  on  what  is  the  normal  synoptic  sounding.  The  synoptic 
tin  es  are  0000  and  1200  (IM  F  but  balloons  were  released  earlier  than  this,  usually 
Oti.10  and  laviO  li TO.  The  ADAS  sonde  was  hung  about  5  m  below  the  synoptic  sonde 
tnd  thev  were  carried  aloft  bv  1  200  g  balloons  with  a  vertical  ascent  rate  of  approxi¬ 
mately  ">  m  sec.  At  non-synoptic  times,  only  the  ADAS  sonde  was  used,  anil  was 
carried  aloft  by  a  1  00  g  balloon  with  a  vertical  ascent  rate  of  approximately  2  m  sec. 


cl  speed  of  approximately  8  n:  see.  The  profile  was  taken  ■  ."lose  -  po.-.-.ibh 
in  space  and  time  to  the  balloon  releases.  The  time  factor  for  synoptic  :  •  s 

was  son'etinies  larger  than  desired  because  the  balloon  could  not  be  delayed.  A 
total  of  17  runs  was  made  over  a  10  day  period.  The  times  are  given  as  Pastern 
Daylight  Time  (DDT'.  Table  1  lists  the  dates  and  times  of  the  various  runs. 

Table  1.  Schedule  of  Kvents 


No. 

Date 

A  eft 

ADAS 

NWS 

1 

9/  17 

0653-0705 

0653-0705 

2 

9/18 

0659-0727 

0649-0700 

0649-0700 

3 

9/18 

1452-1521 

1510-1529 

4 

9/19 

0643-07 14 

0648-0659 

0648-0659 

5 

9/20 

0701-0735 

0636-0647 

0636-0647 

6 

9/20 

1052- 1107 

0934-0955 

7 

9/21 

0724-0742 

0630-0641 

0630-0641 

3 

9/23 

0632-0644 

0632-0644 

9 

9/24 

0710-0725 

0626-0637 

10 

9/25 

0642-0648 

0636-0645 

0636-0645 

1 1 

9/25 

0932-0955 

0930-0947 

12 

9/25 

1606-1623 

1526-1544 

13 

9/25 

1906-1917 

1906-1917 

14 

9/26 

0631-0648 

0638-0650 

0638-0650 

15 

9/26 

0853-0903 

0811-0830 

16 

9/26 

1407-1423 

1435-1452 

17 

9/26 

1833-1844 

1833-1844 

4.  DATA 

Doth  radiosondes  measure  pressure  (P),  temperature  (T),  and  relative 
humidity  (KM)  as  commutated  functions  of  time.  As  mentioned  earlier,  tile  rate 
of  commutation  for  the  two  sondes  is  different.  The  height  to  which  these  measure¬ 
ments  are  assigned  is  computed  using  the  surface  height  and  the  hydrostatic  equation 

dP  =  pgdz  (2) 


A  A  .>  A  AAA  .  -  A  l'  «  -  A.VA  A  .V.VA  A  A  A  A  .  -  .  ■  A  A  A  A  .  •  A  -  A  A  A  .  *  A  .  A  A  A  .  -  A  A  A  .  *  A  A  A  -  'A  AA 
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where 


p  is  air  density. 


g  is  acceleration  due  to  gravity,  and 


z  is  altitude. 


The  index  of  refraction  is  computed  from 
48 10  e 


N  -  ^ 


-) 


where 


Cl) 


T  is  temperature  in  degrees  Kelvin, 

1’  is  total  atmosphere  pressure  in  mb, 
e  is  partial  pressure  of  water  vapor  in  mb, 
and  e  is  calculated  using  the  empirical  relationship 

ew  =  6-1073  '  IOe*e  ( u\.V±r)  ■ 


The  modified  refractive  index  is  then  calculated  from 


(4) 


M  =  N  +  0.  157  Z 


where  Z  is  m  meters. 

In  performing  these  calculations,  all  of  the  data  from  the  ADAS  sonde  were 
iv.* liable  to  us  while  only  the  "significant"  data  from  the  synoptic  soundings  were 
used.  I'his  was  partly  by  choice  since  in  our  field  experiments  on  propagation 
paths,  only  these  data  would  be  available. 

figure  3  shows  the  air  temperature  and  dew  point  temperature  from  a  repre¬ 
sentative  sounding  taken  on  26  September  at  0630  I.DST.  Dew  point  temperature  is 
calculated  from 


T 


D 


24  1 .  0  log  -  (e  id.  107  8) 
_ "  1  0  w _ 

7.63  -  log  (e  ,6.  107 8) 


(til 


0638  EDT  26  SEP  85 


,HRTHRM  MR 


TEMPERATURE  (°C> 


Figure  3.  Atmospheric  nemperature  (T)  and  Dew  Point  Temperature  (Tp>t 

r  • _  r\ _ i  f).,  _ _ 1  _  m:-  lj.  rtr»oo  or?  o _ * _ i _ .  _  i  n  r 


From  Dual  Radiosonde  Flight  0638,  26  September  1985 


Figure  4  shows  the  N  and  M  piofiles  for  this  same  sounding.  Since  the  refiactom- 
eter  reads  out  directly  in  refractive  index,  all  that  is  required  is  the  conversion 
of  the  reading  to  M  units,  using  Eq.  (51.  Figure  5  shows  the  aircraft  sounding 
taken  closest  in  time  to  the  balloon  sounding  of  Figure  4.  The  balloon  sounding 
data  are  included  in  the  figure  and  the  profiles  are  offset  from  each  other  bv  40  M 
units  to  avoid  confusion.  A  radio  duct  is  indicated  in  this  figure,  it  is  quite 
apparent  that  the  modified  refractive  index  (Ml  is  a  much  better  indicator  of  ducting 
than  \.  Also  shown  by  way  of  definition  are  the  duct  thickness,  the  ducting  layer 
(dW/dZ  <  0)  and  the  AM. 


HEIGHT  (m)  HEIGHT 


0G38  EDT  26  SEP  85 


CHRTHRM  MR 


DUCT 

DUCTING 
LAYER  ! 

HICKNESS 

*  i 

H  AM  k/ 

V  /  / 


- RCF 

- RDf* 

- NWS 


I 


100  M  units 


MODIFIED  REFRRCTIVE  INDEX  M 

Figure  5.  Modified  Refractive  Index  (M)  Profiles  for  Refractometer  and 
Balloon  Data  of  0638,  26  September  1985 
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MODIFIED  REFRRCTIVE  INDEX  M 

Figure  6.  Modified  Refractive  Index  (M)  Profiles  for  Events  1  Through  6 
September  1985 
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5.  DISCUSSION 


In  comparing  the  two  types  of  radiosonde  soundings,  it  was  our  intent  to  under¬ 
stand  the  products  of  the  two  as  we  might  use  them.  Our  interest  is  in  the  refrac¬ 
tive  index  profiles  for  use  in  studying  propagation  paths.  When  we  use  data  from 
the  ADAS  system,  we  have  access  to  all  the  data.  When  we  use  data  from  the 
synoptic  NWS  radiosonde  we  are  privy  to  only  the  transmitted  "significant  level" 
data.  The  results  should  not  be  construed  as  implying  anything  about  the  synoptic- 
usefulness  of  either  instrument. 

In  general,  the  agreement  between  the  two  radiosonde  data  bases  is  excellent. 
When  the  two  sets  of  atmospheric  temperatures  are  compared  at  the  significant 
levels  as  reported  from  the  synoptic  sonde,  the  mean  difference  is  found  to  be 
0.  033°  C  and  the  standard  deviation  of  the  differences,  0.  44°C.  When  the  sets  of 
dew  point  temperatures  are  similarly  compared,  the  mean  difference  is  found  to 
be  0.  123°  C  and  the  standard  deviation  of  the  differences,  1.  85°C.  In  this  second 
comparison,  data  were  omitted  when  both  dew  point  depressions  were  greater  than 
20°  C  or  when  the  synoptic  radiosonde  reported  a  30°  C  depression.  Dew  point 
temperature  difference  in  these  dry  conditions  is  not  important  to  refractive  index. 
One  would  expect  good  agreement  since  both  sounders  use  the  same  humidity  sensor, 
a  carbon  humidity  element  (VIZ  Mfg.  No.  1286-  162),  and  the  same  temperature 
measuring  technique,  a  thermistor  with  a  white  paint  coating  exposed  directly  to 
the  air.  The  agreement  is  visible  in  Figure  3  where  temperature  and  dew  point 
temperature  are  plotted  as  functions  of  altitude.  In  addition  to  the  vertical  profiles 
of  temperature  from  the  radiosondes,  a  vertical  profile  of  temperature  from  the 
aircraft  was  also  available  and  can  be  seen  to  agree  with  the  radiosonde  data  to 
better  than  0.  5°C  (Figure  9).  An  aerodynamic  heating  correction  of  1.0°  C  was 
made  to  the  aircraft  measurement. 

Since  we  use  the  data  to  determine  atmospheric  effects  on  microwave  communi¬ 
cation  systems,  it  is  important  to  compare  the  refractive  index  and  refractive  index 
gradients  as  measured  or  derived  from  these  systems.  To  do  this  we  will  use  the 
modified  refractive  index  (Ml  profiles  since  this  allows  easy  visual  identification  of 
radio  ducting  situations,  dM/dZ  <  0.  The  plots  of  MvsZ  for  cases  including  aircraft 
flights  show  that  aircraft  data  indicate  the  greatest  number  of  ducting  situations 
and  the  NWS  sonde  the  least  number.  If  we  ignore  the  first  50  m  above  the  sur¬ 
face,  there  is  only  one  instance  where  a  duct  thicker  than  50  m  is  indicated  in  the 
ADAS  sonde  data  and  not  in  the  aircraft  data.  This  appears  in  the  18  September 
1452  flight,  as  shown  in  Figure  6  as  event  No.  3,  at  an  altitude  of  600  m.  The 
duct  is  only  80  m  thick  and  the  soundings  are  18  min  apart.  Therefore  the  dis¬ 
crepancy  is  probably  atmospheric  and  not  instrumental  in  origin. 
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Figure  9.  Atmospheric  Temperature  Profiles  for  Aircraft  and 
Radiosondes  for  19  September  1985  Event  No.  2 


In  examining  the  events  where  the  ref ractometer  and  the  ADAS  sonde  are  less 
than  20  min  apart  we  found  four  instances  where  the  aircraft  data  showed  a  duct 
greater  than  50  m  thick  that  was  not  indicated  by  the  radiosonde.  Two  additional 
ducts,  although  indicated  by  the  ADAS  sonde,  were  indicated  as  less  than  50  m 
thick.  In  no  case  was  a  duct  thicker  than  60  m  indicated  by  the  aircraft  and  not 
indicated  by  the  ADAS.  Of  the  11  ducts  indicated  by  both  the  ref  ractometer  and  the 
ADAS,  the  refractometer  shows  a  larger  AM  in  8  cases  and  the  ADAS  AM  is  larger 
in  2  cases. 

The  M  profile  data  from  the  NWS  sonde  does  not  show  ducting  in  several  in¬ 
stances  where  both  the  aircraft  and  ADAS  data  indicate  the  ducting  criterion  has 
been  satisfied.  Two  of  the  most  notable  occurrences  are  in  event  No.  14,  shown  in 
Figure  8,  and  event  No.  2,  shown  in  Figure  6,  at  altitudes  of  800  m  and  500  m, 
respectively.  When  there  are  just  balloon  releases,  there  are  additional  instances 
where  ducting  is  indicated  by  the  ADAS  data  and  not  bv  the  NWS  data.  One  of  the 
most  dramatic  cases  was  event  No.  17  (Figure  8)  on  26  September  1985  at  1888. 

The  ADAS  instrument  shows  a  strong  duct  between  600  and  700  m,  a  second  strong 


duct  between  800  and  1050  m,  and  a  third  weaker  duct  near  1300  m.  The  synoptic 
data  show  only  the  first  of  these.  The  temperatures  and  dew  point  temperatures 
for  this  flight,  illustrated  in  Figure  10,  show  that  the  synoptic  data  do  not  indicate 
the  two  tongues  of  dry  air  that  caused  the  ducting.  The  layer  at  1000  m,  at  least, 
was  missed  because  of  the  way  the  data  are  commutated  in  the  NWS  sonde.  The 
dry  layer  occurred  where  the  reference  contact  caused  the  humidity  contacts  to 
be  about  60  sec  apart. 
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TEMPERRTURE  (°C) 

Figure  10.  Atmospheric  Temperature  (T)  and  Dew  Point  Temperature  (T.J 
for  Balloon  Sounding  of  1833,  26  September  1985 

Examination  of  all  the  events  involving  both  the  ADAS  and  the  synoptic  sonde 
reveals  11  instances  where  the  ADAS  showed  a  duct  greater  than  50  m  thick  and 
the  synoptic  data  did  not  indicate  ducting.  Of  these,  7  were  greater  than  100  m 
and  only  1  was  greater  than  150  m.  In  the  10  instances  where  both  the  ADAS  and 
the  synoptic  sonde  show  ducting,  the  ADAS  sonde  indicated  a  larger  value  of  AM 


6.  CONCLUSIONS 


The  higher  commutation  rate  of  the  ADAS  instrument  makes  it  possible  to 
detect  all  ducts  more  than  60  m  thick  when  used  on  balloons  with  ascent  rates  of 
6  m/sec  in  the  first  few  kilometers.  In  one  case,  the  data  from  the  synoptic  sonde 
failed  to  indicate  a  duct  205  m  thick,  (event  No.  17).  In  general,  the  synoptic  data 
will  detect  all  ducts  more  than  250  m  thick  and  only  rarely  miss  ducts  more  than 
150  m  thick.  When  all  data  are  available  from  an  instrument  with  a  higher  commu 
tation  rate  (5  sec  frame  rate),  any  duct  more  than  60  m  thick  should  be  detected. 

In  our  propagation  experiments  we  have  used  ADAS  instruments  and  rise  rates  of 
2.  5  m/sec,  so  that  ducts  as  thin  as  30  m  are  detectable. 

The  aircraft  data  is  limited  by  the  update  rate  of  the  refractive  channel,  1  sec 
and  the  rise  rate  of  the  aircraft,  in  our  tests  about  3  m/sec.  This  results  in  reso¬ 
lution  of  about  6  m  ducts. 

In  addition  to  missing  some  of  the  thinner  ducts,  both  radiosonde  instruments 
underestimate  the  change  in  M  in  the  ducts  relative  to  the  aircraft  refractometer. 
This  has  two  types  of  effects;  the  duct  thickness  is  less,  and  the  critical  angle 
necessary  for  a  radio  wave  to  escape  is  greater.  There  is  no  indication  of  either 
instrument  indicating  ducting  erroneously. 

These  tests  would  indicate  that  any  study  of  the  frequency  of  occurrence  and 
strength  of  radio  ducting  based  on  synoptic  data  will  underestimate  the  frequency 
for  ducts  with  thicknesses  less  than  300  m  and  will  underestimate  the  strength  of 
the  ducts  in  general.  The  limited  amount  of  our  testing  is  not  sufficient  to  make 
quantitative  estimates  of  this  deficiency. 


